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In  this  study,  a  large-scale  integrated  modeling  system  (IMS)  was  developed  for  supporting  climate- 
change  impact  analysis  and  adaptation  planning  under  multi-level  uncertainties.  A  number  of 
methodologies  were  seamlessly  incorporated  within  IMS,  including  fuzzy-interval  inference  method 
(FIIM),  inexact  energy  model  (IEM),  and  uncertainty  analysis.  The  system  could  (i)  encompass  multi¬ 
ple  technologies,  energy  resources,  and  sub-sectors,  and  climate  change  impact  analysis  into  a  general 
modeling  framework,  (ii)  address  interactions  of  climate  change  impacts  on  multiple  energy  sub-sectors 
and  resources  within  an  EMS,  (iii)  identify  optimal  adaptation  strategies  of  an  EMS  to  climate  change 
impact  through  a  two-step  procedure,  (iv)  deal  with  multiple  levels  of  uncertain  information  associated 
with  processes  of  climate  change  impact  analysis  and  adaptation  planning,  and  (v)  seamlessly  combine 
climate  change  impact  analysis  results  with  inexact  adaptation  planning.  It  could  provide  decision  mak¬ 
ers  a  comprehensive  view  on  the  EMS  as  well  as  the  corresponding  adaptation  schemes  under  climate 
change,  greatly  improving  the  robustness  and  completeness  of  the  decision-making  processes.  The  gen¬ 
erated  solutions  can  provide  desired  energy  resource/service  allocation  with  a  minimized  system  cost, 
a  maximized  system  reliability  and  a  maximized  energy  security  under  varied  climate  change  impact 
levels,  as  well  as  multiple  levels  of  uncertainties.  In  a  companion  paper,  the  developed  method  is  applied 
to  a  real  case  for  the  long-term  planning  of  waste  management  in  the  Province  of  Manitoba,  Canada. 
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1.  Introduction 

Sustainable  energy  utilization  has  become  a  prominent  policy 
target  in  many  countries  across  the  world  [1,2].  However,  there  are 
many  concerns  over  the  adoption  of  various  energy  resources,  such 
as  energy-shortage  risk,  greenhouse  gas  emission,  and  renewable 
energy  development,  posing  great  challenges  for  decision  mak¬ 
ers  in  multi-scale  energy  management  systems  (EMSs).  Also,  the 
management  efforts  are  complicated  by  climate  change  due  to 
its  significant  effects  on  statistics  of  many  energy  sources  (e.g., 
max,  min,  mean  values  and  variance  of  hydropower,  solar  and 
wind  energies).  Moreover,  such  efforts  are  subject  to  multiple 
levels  of  uncertainties  that  are  associated  with  the  forecast  of 
changing  climatic  conditions  (e.g.,  uncertain  shifts  of  precipitation 
patterns),  the  interactions  between  energy  resources/technologies 
and  climatic  variations  (e.g.,  uncertain  fluctuation  of  hydropower 
generation  in  responding  to  precipitation  variations),  and  the  plan¬ 
ning  of  EMSs  (e.g.,  uncertain  estimation  of  energy  prices).  This  leads 
to  a  variety  of  complexities  in  relevant  decision-making  under  cli¬ 
mate  change.  These  complexities  may  be  further  multiplied  by 
multi-period,  multi-facility  and  multi-objective  features  as  well  as 
system  dynamics.  Thus,  system  analysis  techniques  are  desired  to 
assist  in  formulating  long-term  energy  management  and  allocation 
plans,  which  will  be  helpful  not  only  for  analyzing  tradeoffs  among 
various  socio-economic,  energy-related,  and  environmental  objec¬ 
tives,  but  also  for  providing  cost-effective  adaptation  schemes  in 
responding  to  varied  impact  levels  of  climate  change. 

Previously,  a  plenty  of  efforts  were  undertaken  for  investigat¬ 
ing  impacts  of  climate  change  on  multi-scale  energy  management 
systems  and  planning  the  corresponding  adaptation  capacities 
[3-13,74-78].  For  instance,  Peramunetilleke  [5]  assessed  impacts 
of  climate  change  on  Sri  Lanka’s  hydroelectric  production.  They 
believed  that  power  supply  from  hydropower  without  back-up 
sources  would  be  highly  risky  under  extreme  climatic  events. 
Allen  and  Christensen  [14]  briefly  reviewed  the  impacts  of  cli¬ 
mate  change  and  suggested  that  a  diverse  range  of  issues  would 
contribute  to  this  process.  Sutherland  [15]  qualitatively  analyzed 
potential  impacts  of  climate  change  on  a  number  of  energy- 
intensive  industries  in  US.  Breslow  and  Sailor  [16]  investigated 
impacts  of  climate  change  on  wind  speeds  and,  hence,  on  wind 
power  generation  in  US.  A  series  of  GCM  outputs  from  the  Cana¬ 
dian  Climate  Center  and  the  Hadley  Center  were  used  for  providing 
a  range  of  possible  variations  in  seasonal  mean  wind  magnitude. 
Potential  impacts  of  climate  change  on  heating  and  cooling  energy 
demands  were  investigated  through  building  energy  simulation 
models  under  hourly  weather  scenarios  and  were  applied  to  the 
Zurich-Kloten  area  of  Switzerland  [17].  Gaterell  and  McEvoy  [10] 
examined  uncertain  effects  of  climate  change  on  performance  of 
energy  insulation  measures.  Christenson  et  al.  [8]  investigated 
impacts  of  climate  change  on  energy  demand  of  residential  build¬ 
ings  through  the  degree-days  method.  A  procedure  to  estimate 
heating  and  cooling  degree-days  from  monthly  temperature  data 
was  proposed,  verified,  and  applied  to  four  representative  buildings 
in  Switzerland.  Mirasgedis  et  al.  [18]  adopted  a  regional  climate 
model  (i.e.,  PRECIS)  to  predict  future  climatic  conditions  under  dif¬ 
ferent  emissions  scenarios.  The  predicted  data  were  used  as  inputs 
to  a  multiple-regression  model  to  examine  sensitivity  of  electricity 
demand  in  Greece’s  interconnected  power  system  under  various 
climatic  and  socio-economic  conditions.  Pryor  and  Barthelmie  [19] 
qualitatively  assessed  impacts  of  climate  change  on  wind  energy 
sources,  particularly  on  operation  and  maintenance  of  several  wind 
farms.  Based  on  a  qualitative  method,  Whitmarsh  [20]  investigated 
various  responses  to  the  impacts  of  climate  change  in  many  energy 
management  systems.  Guan  [21]  attempted  to  evaluate  impacts 
of  climate  change  on  building  environment  through  combination 
of  weather-data  forecasting  and  building  simulation  techniques. 


Bassi  and  Baer  [22]  employed  an  integrated  modeling  approach 
to  carry  out  a  country-wide,  cross-sector  analysis  of  the  interac¬ 
tions  among  Ecuador’s  energy,  social,  economic  and  environmental 
sectors  under  climate  change.  Isaac  and  van  Vuuren  [23]  assessed 
potential  energy  consumptions  for  future  residential  heating  and 
air  conditioning  under  climate  change. 

Under  a  series  of  climate  change  impact  levels,  adaptation  plans 
of  EMSs  need  to  be  systematically  considered  by  decision  mak¬ 
ers.  Over  the  past  two  decades,  a  large  number  of  studies  were 
undertaken  on  adaptation  planning  toward  climate  change  in  envi¬ 
ronmental  and  energy  management  systems  planning  [24-50]. 
For  example,  Huang  et  al.  [51]  proposed  an  inexact-fuzzy  multi¬ 
objective  programming  model  for  supporting  adaptation  planning 
of  land  resources  management  in  Mackenzie  Basin  under  changing 
climatic  conditions.  Stakhiv  [52]  suggested  that  policy  makers  and 
water  resource  managers  should  be  aware  of  the  evolving  informa¬ 
tion  on  climate  change  impacts  as  an  activity  that  was  preparatory 
to  sound  decision  making  on  current  water  resource  management 
actions.  Tol  et  al.  [24]  made  a  thorough  investigation  on  studies 
related  to  climate  change  adaptation.  They  suggested  that  there 
should  be  four  categories  of  adaptation  measures  toward  climate 
change,  including:  no  adaptation,  arbitrary  adaptation,  observed 
adaptation,  and  modeled  adaptation  (optimization).  Since  climate 
change  poses  significant  challenges  for  water  resource  manage¬ 
ment  in  Canada,  de  Loe  et  al.  [53]  discussed  issues  related  to 
identification  of  adaptation  measures  to  water  supply  systems 
in  the  near  future.  Naess  et  al.  [33]  examined  the  adaptation  of 
infrastructural  institutions  under  climate  change  in  Norway.  Two 
municipalities  were  used  as  examples  for  investigating  institu¬ 
tional  responses  to  floods  under  climate  change.  Dessai  and  Hulme 
[54]  presented  an  assessment  framework  for  identifying  adap¬ 
tation  strategies  that  were  robust  (i.e.,  insensitive)  to  climate 
change  uncertainties.  The  framework  was  applied  to  a  real-world 
case  for  supporting  water  resources  management  in  eastern  Eng¬ 
land.  Osbahr  et  al.  [55]  explored  cross-scale  dynamics  in  coping 
and  adjusting  responses  toward  climate  change  based  on  qualita¬ 
tive  data  from  a  real-world  case  in  Mozambique.  More  recently, 
Hoffmann  et  al.  [56]  empirically  examined  the  determinants  of 
adaptation  measures/policies  toward  climate  change.  They  also 
analyzed  the  influence  of  possible  climate  change  effects  on  the 
scope  of  adaptation  measures  through  a  number  of  linear  regres¬ 
sion  models.  Bryan  et  al.  [41  ]  analyzed  the  adaptation  strategies  of 
farmers  in  several  countries  and  analyzed  the  factors  influencing 
their  decisions.  Deressa  et  al.  [57]  identified  major  strategies  used 
by  farmers  to  adapt  to  climate  change  in  Nile  Basin  of  Ethiopia.  They 
also  analyzed  the  factors  that  might  affect  their  choices  of  method 
and  the  barriers  to  adaptation. 

Particularly,  a  number  of  studies  were  undertaken  for  exam¬ 
ining  potential  adaptation  strategies  of  energy  sectors  in  many 
countries  under  climate  change.  For  instance,  Simeonova  [58]  com¬ 
plied  and  analyzed  a  variety  of  policies  and  measures  for  adapting 
to  climate  change  in  the  energy  sectors  of  many  central  and  east¬ 
ern  European  countries.  Matondo  et  al.  [59]  examined  the  impacts 
of  climate  change  on  hydrological  regimes,  water  resources  man¬ 
agement,  and  hydropower  generation.  Ruth  and  Lin  [60]  explored 
potential  impacts  of  climate  change  on  consumptions  of  natural 
gas,  electricity,  and  heating  oil  by  residential  and  commercial  sec¬ 
tors  in  the  state  of  Maryland  of  the  US.  Time  series  analysis  was 
then  used  to  quantify  relationships  between  historical  tempera¬ 
ture  and  energy  demand.  A  dynamic  computer  model  based  on 
those  relationships  was  used  for  simulating  future  energy  demand 
under  a  range  of  energy  prices,  temperatures,  and  other  driving 
factors.  Mansur  et  al.  [61  ]  adopted  a  national  energy  model  of  fuel 
choice  to  generate  climate  change  adaptation  strategies  at  both 
households  and  firm  scales.  They  suggested  that  increasing  elec¬ 
tricity  consumption  for  space  cooling  would  be  observed  with  a 


Y.P.  Cai  et  al.  /  Renewable  and  Sustainable  Energy  Reviews  15  (201 1 )  2779-2790 


2781 


reduction  in  energy  consumption  for  space  heating  under  climate 
change.  Overall,  energy-related  expenditures  would  likely  increase 
in  America  due  to  the  impacts  of  climate  change  on  residential 
heating  and  cooling,  resulting  in  potential  welfare  damages.  Jenk¬ 
ins  et  al.  [62]  attempted  to  quantify  how  climate  change  will  have  a 
direct  effect  on  heating  and  cooling  energy  use  in  future  office  envi¬ 
ronments.  The  results  confirmed  the  importance  of  demand-side 
management  before  assessing  the  supply-side  opportunities  under 
extreme  climatic  events.  The  study  also  highlighted  the  importance 
and  possibilities,  of  adaptation  strategies  to  future  climates,  and 
the  benefits  of  promoting  heating-dominated  buildings  instead  of 
cooling-dominated  ones.  Weisser  et  al.  [63]  used  a  qualitative  anal¬ 
ysis  method  to  explore  the  conditions  under  which  nuclear  power 
could  adapt  to  climate  change  over  a  long  term  period. 

Previous  studies  on  climate  change  impact  analysis  related  to 
EMSs  are  mostly  based  on  qualitative  methods,  and,  thus,  could 
merely  generate  qualitative  information.  Many  studies  were  con¬ 
ducted  to  analyze  the  impacts  of  climate  change  on  energy  sectors 
as  well  as  to  examine  the  corresponding  adaptation  strategies 
through  integrating  many  quantitative  methods  and/or  GCM/RCM 
results.  However,  direct  incorporation  of  GCM/RCM  results  into 
climate-change  impact  analysis  for  multi-scale  EMSs  was  limited  to 
one  or  several  climatic  sensitive  parameters  such  as  temperature, 
precipitation,  and  humidity,  which  might  be  unreasonable  in  real¬ 
istic  problems.  In  addition,  most  of  the  previous  studies  focused  on 
impacts  of  climate  change  on  an  individual  energy  sub-sector,  pro¬ 
cess  or  technology  (e.g.,  residential,  commercial  and  hydropower 
generation  sub-sectors).  There  was  a  lack  of  studies  that  could 
comprehensively  reflect  and  analyze  impacts  of  climate  change  on 
multiple  energy  sub-sectors  and  technologies,  and  thus  analyze  the 
impacts  from  an  entire  EMS  point  of  view.  Moreover,  the  process 
of  impact  analysis  involved  multiple  levels  of  imprecise  informa¬ 
tion,  which  was  largely  neglected  in  the  previous  studies.  In  fact, 
the  process  of  climate  change  impact  analysis  is  subject  to  a  vari¬ 
ety  of  uncertainties.  For  example,  in  most  of  the  quantitative  or 
semi-quantitative  methods  based  on  GCM/RCM  results,  there  are 
multiple  formats  of  uncertainties  embedded  within  not  only  the 
input  information,  but  also  the  weather-forecasting  mechanisms. 
Furthermore,  impact  analysis  could  merely  provide  answers  to  the 
question  of  “which  sector/process  within  an  energy  management 
system  would  be  the  most  vulnerable  under  climate  change”.  Deci¬ 
sion  makers  may  need  to  know  “if  there  is  any  necessity  to  take 
adaptation  actions”  and  “what  the  optimal  adaptation  actions  for 
the  energy  sector  under  multiple  system  uncertainties  will  be”.  The 
questions  need  to  be  answered  through  EMSs  adaptation  planning. 
Likewise,  in  the  process  of  adaptation  planning  based  on  optimiza¬ 
tion  models,  relevant  parameters  and  coefficients  may  be  subject 
to  multiple  formats  of  uncertainties  (e.g.,  intervals,  probabilistic 
and  possibilistic  distributions).  Until  now,  there  were  no  reported 
studies  that  could  link  adaptation  planning  with  climate-change 
impact  analysis  and  could  address  multi-level  uncertainties  asso¬ 
ciated  with  impact  analysis  as  well  as  the  subsequent  adaptation 
planning  process.  Therefore,  an  integrated  modeling  system  for 
seamlessly  coupling  climate  change  impacts  with  EMSs  adapta¬ 
tion  planning  is  desired,  particularly  for  a  region  that  is  heavily 
dependent  on  renewable  energy  resources  and  is,  thus,  vulnerable 
to  climate  change. 

2.  Methodology 

2.2.  A  systematic  perspective  on  adaptation  planning  under 
climate  change 

Almost  all  energy-related  processes  and  factors  would  be 
directly  and/or  indirectly  influenced  by  changing  climatic  condi- 


Fig.  1.  Climate  change  impacts  and  adaptation  schemes  of  EMSs. 


tions  (Fig.  1).  Among  them,  end-user  demands,  renewable  energy 
resources  and  the  corresponding  utilization  facilities  would  prob¬ 
ably  be  the  most  sensitive  and  vulnerable  ones.  For  instance, 
availabilities  of  renewable  energy  resources  (e.g.,  hydropower  and 
wind  energy)  are  highly  dependent  on  the  statistics  (e.g.,  min, 
peak  and  mean  values  as  well  as  variance)  of  many  meteorolog¬ 
ical  variables  that  might  fluctuate  within  a  certain  range  due  to 
climate  change.  Such  variations  of  renewable  energy  availabilities 
would  then  affect  operating  statuses  of  relevant  facilities,  result¬ 
ing  in  changes  in  their  energy  outputs.  Similarly,  in  responding 
to  variations  of  climatic  factors/conditions,  such  as  temperature, 
solar  radiation,  and  humidity,  end-user  demands  for  space  heat¬ 
ing,  water  heating,  and  space  cooling  would  correspondingly  vary 
over  a  certain  range.  Apparently,  these  impacts  would  not  affect 
an  individual  energy  sub-sector  (e.g.,  energy  demand  and  produc¬ 
tion)  in  a  separate  way.  Instead,  they  interact  with  each  other  and, 
thus,  would  cause  an  integrated  impact  on  the  entire  EMS.  Such 
an  integrated  impact  over  a  long-term  planning  period  needs  to 
be  systematically  analyzed  by  decision  makers  in  order  to  discern 
the  most  optimal  adaptation  schemes  for  the  entire  EMS  under 
climate  change.  However,  such  a  task  is  difficult  due  to  a  vari¬ 
ety  of  uncertainties  associated  with  changing  climatic  conditions 
as  well  as  decision  makers’  subjective  judgments  on  impacts.  For 
example,  due  to  the  lack  of  sufficient  information  related  to  the 
evaluation  guidelines  (e.g.,  impact  levels),  fuzziness  is  always  asso¬ 
ciated  with  the  process  of  impact  identification  and  integration, 
which  is  normally  caused  by  vague  and  imprecise  linguistic  terms. 
Correspondingly,  the  evaluation  of  occurrence  chances  of  climate 
change  and  their  consequences  on  EMSs  is  also  subject  to  uncertain¬ 
ties  and  is  difficult  to  be  handled  through  conventional  methods. 
This  implies  that  planning  activities  for  EMSs  would  be  misled  into 
a  deviated  or  even  false  direction  without  a  comprehensive  consid¬ 
eration  over  the  integrated  impacts  of  climate  change. 

Moreover,  there  is  a  variety  of  uncertain  information  associated 
with  decision  process,  which  can  only  be  expressed  as  interval 
numbers  without  known  distributions.  Such  uncertain  informa¬ 
tion  need  to  be  effectively  incorporated  into  the  planning  process 
and  the  obtained  decision  alternatives,  which  represents  a  differ¬ 
ent  level  of  uncertainties  compared  with  those  implied  in  climate 
change  analysis.  Therefore,  in  this  research,  the  integrated  climate 
change  impacts  will  be  incorporated  into  the  modeling  framework 
to  facilitate  identifying  adaptation  plans  with  increased  robustness 
under  various  changing  conditions.  Specifically,  a  fuzzy-interval 
inference  method  (FIIM)  and  an  inexact  energy  model  (IEM)  will 
be  developed  and  then  integrated  into  a  general  modeling  system 
(IEMS)  as  well  as  applied  to  the  province  of  Manitoba  (IEMS- 
Manitoba,  presented  in  a  companion  paper). 
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Fig.  2.  Membership  function  of  climate  change  impacts  on  energy  supply  (i.e., 
renewable  availability  factors  based  on  renewable  energy  resources).  Note : 
AF  =  availability  factor. 

2.2.  Coupled  climate -change  impact  analysis  under  uncertainty 
2.2.1.  Impacts  on  energy  supply 

Considering  an  EMS  wherein  electricity  production  mainly 
occurs  based  on  renewable  energy  resources  such  as  hydropower 
and  wind  energy,  over  a  long-term  period,  climate  change  will  have 
direct  impacts  on  the  availability  of  these  resources  due  to  potential 
variations  in  many  climatic  factors  and  conditions  such  as  precipi¬ 
tation  types,  radiation  intensities,  and  wind  speeds.  Consequently, 
the  changes  in  resource  availability  would  affect  operation  sta¬ 
tuses  of  relevant  electricity  generation  facilities.  In  this  research, 
an  availability  factor  is  introduced,  which  is  defined  as  the  amount 
of  time  that  a  hydroelectric  power  station  (or  wind  farm)  needs 
to  produce  electricity  over  a  certain  period,  divided  by  the  total 
time  of  this  period.  Then,  five  fuzzy  sets  are  assigned  to  repre¬ 
sent  different  levels  of  impacts  of  climate  change  on  the  availability 
factors  of  renewable-energy-based  power  generation  facilities,  cor¬ 
responding  to  “low  impact”,  “low-to-medium  impact”,  “medium 
impact”,  “medium-to-high  impact”  and  “high  impact”.  The  con¬ 
struction  of  membership  functions  for  these  five  fuzzy  sets  relies 
on  the  expertise  and  experience  of  experts  and  stakeholders.  The 
involved  experts  and  stakeholders  should  have  in-depth  knowl¬ 
edge  about  the  province  and  the  energy  management  system,  and 
can,  thus,  provide  valuable  inputs  for  quantifying  the  uncertainties 
[64,65].  Fuzzy  logic  is  then  employed  for  processing  such  subjective 
opinions. 

In  this  study,  fuzzy  sets  corresponding  to  different  levels  of 
climate  change  impacts  were  obtained  through  questionnaire  sur¬ 
veys.  Forty-eight  persons  were  interviewed  about  their  opinions  on 
the  corresponding  relations  between  climate  change  impact  lev¬ 
els  and  renewable  energy  availability  factor  decrements.  As  listed 
in  Table  1,  54.10%  of  the  surveyed  respondents  believed  that  “low 
impact”  should  correspond  to  the  statement  “the  availability  fac¬ 
tor  of  a  renewable-energy-based  power  generation  facility  will 
decrease  by  5.00%  or  less”;  58.40%  of  the  respondents  selected:  “the 
availability  factor  of  a  renewable-energy-based  power  generation 
facility  will  decrease  by  approximately  10%”  as  “low-to-medium 
impact”;  58.40%  of  the  respondents  agreed  that  “the  availabil¬ 
ity  factor  of  a  renewable-energy-based  power  generation  facility 
will  decrease  by  around  15.00%”  represents  “medium  impact”; 
60.40%  of  the  respondents  matched  “the  availability  factor  of  a 
renewable-energy-based  power  generation  facility  will  decrease 
by  approximately  20.00%”  with  “medium-to-high  impact”;  and 
77.10%  of  the  respondents  considered  “the  availability  factor  of  a 
renewable-energy-based  power  generation  facility  will  decrease 
25.00%  or  more”  as  “high  impact”.  According  to  Hwang  and  Chen 
[66],  the  membership  functions  of  these  fuzzy  sets  can  then  be 


Table  1 

Impacts  of  climate  change  on  energy  supply  (power  generation  from  renewable 
energy  resources). 


The  availability  factor  of  renewable 
energy  resources  will  decrease  by 
approximately: 

No.  of  respondents 

Percentage  (%) 

Low  impacts 

1%  or  less 

0 

0.00 

2%  or  less 

0 

0.00 

3%  or  less 

13 

27.10 

4%  or  less 

9 

18.80 

5%  or  less 

26 

54.10 

6%  or  less 

0 

0.00 

7%  or  less 

0 

0.00 

8%  or  less 

0 

0.00 

9%  or  less 

0 

0.00 

Total  no.  of  respondents 

48 

100.00 

Low-to-Medium  impacts 

5% 

0 

0.00 

6% 

0 

0.00 

7% 

0 

0.00 

8% 

10 

20.80 

9% 

10 

20.80 

10% 

28 

58.40 

15% 

0 

0.00 

20% 

0 

0.00 

25% 

0 

0.00 

Total  no.  of  respondents 

48 

100.00 

Medium  impacts 

5% 

0 

0.00 

6% 

0 

0.00 

7% 

0 

0.00 

8% 

0 

0.00 

9% 

10 

20.80 

10% 

10 

20.80 

15% 

28 

58.40 

20% 

0 

0.00 

25% 

0 

0.00 

Total  no.  of  respondents 

48 

100.00 

Medium-to-High  impacts 

5% 

0 

0.00 

6% 

0 

0.00 

7% 

0 

0.00 

8% 

0 

0.00 

9% 

0 

0.00 

10% 

6 

12.50 

15% 

13 

27.10 

20% 

29 

60.40 

25% 

0 

0.00 

Total  no.  of  respondents 

48 

100.00 

High  impacts 

5%  or  greater 

0 

0.00 

6%  or  greater 

0 

0.00 

7%  or  greater 

0 

0.00 

8%  or  greater 

0 

0.00 

9%  or  greater 

0 

0.00 

10%  or  greater 

0 

0.00 

15%  or  greater 

0 

0.00 

20%  or  greater 

11 

22.90 

25%  or  greater 

37 

77.10 

Total  no.  of  respondents 

48 

100.00 

constructed  based  on  above  investigation  results.  The  constructed 
fuzzy  sets  are  displayed  in  Fig.  2.  In  this  figure,  “L”,  “L-M”,  “M”, 
“M-H”  and  “H”  represent  low,  low-to-medium,  medium,  medium- 
to-high,  and  high  impacts,  respectively. 
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Table  2 

Climate  change  impacts  on  end-user  energy  demands  (heating  and  cooling). 


Energy  demands  for  heating  and 
cooling  will  increase  by  approximately: 

No.  of 

respondents 

Percentage 

(%) 

Low  impacts 

1%  or  less 

23 

47.90 

1.5%  or  less 

14 

29.20 

2.0%  or  less 

11 

22.90 

2.5%  or  less 

0 

0.00 

3.0%  or  less 

0 

0.00 

3.5%  or  less 

0 

0.00 

4.0%  or  less 

0 

0.00 

4.5%  or  less 

0 

0.00 

5.0%  or  less 

0 

0.00 

Total  no.  of  respondents 

48 

100.00 

Low-to-Medium  impacts 

1% 

0 

0.00 

1.5% 

6 

12.50 

2.0% 

19 

39.60 

2.5% 

15 

31.20 

3.0% 

7 

14.60 

3.5% 

0 

0.00 

4.0% 

0 

0.00 

4.5% 

0 

0.00 

5.0% 

0 

0.00 

Total  no.  of  respondents 

48 

100.00 

Medium  impacts 

1% 

0 

0.00 

1.5% 

0 

0.00 

2.0% 

2 

4.20 

2.5% 

10 

20.80 

3.0% 

21 

43.80 

3.5% 

11 

22.90 

4.0% 

4 

8.30 

4.5% 

0 

0.00 

5.0% 

0 

0.00 

Total  no.  of  respondents 

48 

100.00% 

Medium-to-High  impacts 

1% 

0 

0.00 

1.5% 

0 

0.00 

2.0% 

0 

0.00 

2.5% 

1 

2.00 

3.0% 

8 

16.70 

3.5% 

14 

29.20 

4.0% 

21 

43.80 

4.5% 

4 

8.30 

5.0% 

0 

0.00 

Total  no.  of  respondents 

48 

100.00 

High  impacts 

1%  or  greater 

0 

0.00 

1.5%  or  greater 

0 

0.00 

2.0%  or  greater 

0 

0.00 

2.5%  or  greater 

0 

0.00 

3.0%  or  greater 

0 

0.00 

3.5%  or  greater 

1 

2.00 

4.0%  or  greater 

3 

6.30 

4.5%  or  greater 

21 

43.80 

5.0%  or  greater 

23 

47.90 

Total  no.  of  respondents 

48 

100.00 

2.2.2.  Impacts  on  energy  demand 

Similar  to  the  availabilities  of  renewable  energy  resources, 
energy  demands  are  also  sensitive  to  the  variations  in  climatic 
conditions,  especially  those  for  heating  and  cooling  in  residential 
and  commercial  sectors  [67,68,73].  Similar  to  the  identification  of 
climate  change  impacts  on  energy  supply,  the  impacts  on  end- 
user  energy  demands  are  also  categorized  into  five  levels  (i.e., 
low,  low-to-medium,  medium,  medium-to-high,  and  high)  which 


Fig.  3.  Membership  function  of  the  climate  change  impacts  on  end-user  energy 
demand. 

are  presented  as  five  fuzzy  sets.  A  survey  was  conducted  to 
obtain  the  associated  membership  functions.  The  energy  demands 
were  assumed  to  increase  by  1.00-5.00%  compared  to  the  orig¬ 
inal  amounts.  The  survey  results,  listed  in  Table  2  suggest  that: 
“energy  demands  for  heating  and  cooling  in  the  residential  and 
commercial  sectors  will  increase  by  approximately  1.00%  or  less” 
corresponds  to  “low  impact”;  “energy  demands  for  heating  and 
cooling  in  residential  and  commercial  sectors  will  increase  by 
approximately  1.50%”  corresponds  to  “low-to-medium  impact”; 
“energy  demands  for  heating  and  cooling  in  residential  and  com¬ 
mercial  sectors  will  increase  by  approximately  2.00%”  corresponds 
to  “medium  impact”;  “energy  demands  for  heating  and  cooling 
in  residential  and  commercial  sectors  will  increase  by  approxi¬ 
mately  2.50%”  corresponds  to  “medium-to-high  impact”;  “energy 
demands  for  heating  and  cooling  in  residential  and  commercial  sec¬ 
tors  would  increase  approximately  3.00%  or  more”  corresponds  to 
“high  impact”.  Fuzzy  membership  functions  associated  with  these 
five  impact  levels  are  established  according  to  the  above  investiga¬ 
tion  results  (Fig.  3). 

2.2.3.  Coupled  impacts  on  energy  management  systems 

In  order  to  comprehensively  investigate  the  impacts  of  climate 
change  on  the  entire  EMS,  integrated  impact  levels  are  derived 
from  combining  those  impacts  on  both  energy  supply  and  end- 
user  demand  sub-sectors.  However,  it  is  very  hard  to  link  climate 
impacts  on  both  sides  through  conventional  quantitative  methods. 
Therefore,  fuzzy  logic  could  be  an  effective  tool  for  facilitating  the 
quantification  of  such  integrated  impacts,  leading  to  the  develop¬ 
ment  of  a  fuzzy-interval  inference  method  (FIIM). 

In  this  study,  the  determination  of  the  integrated  impact  levels 
will  be  based  on  a  series  of  fuzzy  rules  from  experts  and  stakehold¬ 
ers.  The  impact  levels  were  set  to  include  six  categories,  namely, 
“low”,  “low-to-medium”,  “medium”,  “medium-high”,  “high”  and 
“very-high”.  The  six  levels  are  expressed  as  six  fuzzy  sets.  Then,  the 
fuzzy  logic  operator  “AND”  is  used  to  join  factors  in  the  antecedents 
of  the  rules.  Since  both  environmental  and  health  risks  include 
five  categories  of  fuzzy  events,  there  will  be  a  total  of  150  rules 
(5  x  5  x  6).  Assume  that,  if  a  rule  achieves  the  highest  frequency 
in  the  survey,  then  it  will  be  kept  in  the  rule  base  for  the  deter¬ 
mination  of  the  integrated  impact  level.  It  was  found  that  68.30% 
of  the  surveyed  respondents  selected  “if  both  impacts  on  demand 
and  supply  are  low,  then  the  integrated  impact  will  be  low”;  78.50% 
selected  “if  the  impacts  on  supply  side  is  low  and  demand  side  is 
low-to-medium,  then  the  integrated  impact  level  will  be  low-to- 
medium”;  76.80%  chose  “if  the  impact  on  energy  supply  is  low  and 
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Table  3 

Rules  for  the  coupled  climate  change  impact  level. 


No. 

Antecedent 

Then  the  integrated  impact  is 

If  impacts  on  renewable  energy  availability  factor  is 

And  impacts  on  end-user  demands  is 

L  L-M  M 

M-H 

H 

VH 

1 

Low 

Low 

V 

2 

Low 

Low-to-Medium 

V 

3 

Low 

Medium 

V 

4 

Low 

Medium-to-High 

V 

5 

Low 

High 

V 

6 

Low-to-Medium 

Low 

V 

7 

Low-to-Medium 

Low-to-Medium 

V 

8 

Low-to-Medium 

Medium 

V 

9 

Low-to-Medium 

Medium-to-High 

V 

10 

Low-to-Medium 

High 

V 

11 

Medium 

Low 

V 

12 

Medium 

Low-to-Medium 

V 

13 

Medium 

Medium 

V 

14 

Medium 

Medium-to-High 

V 

15 

Medium 

High 

V 

16 

Medium-to-High 

Low 

V 

17 

Medium-to-High 

Low-to-Medium 

V 

18 

Medium-to-High 

Medium 

V 

19 

Medium-to-High 

Medium-to-High 

V 

20 

Medium-to-High 

High 

V 

21 

High 

Low 

V 

22 

High 

Low-to-Medium 

V 

23 

High 

Medium 

V 

24 

High 

Medium-to-High 

V 

25 

High 

High 

V 

Note :  L  = 

Low;  L-M  =  Low-to-Medium;  M  =  Medium;  M-to-H  =  Medium-to-High;  H  =  High;  VH  =  Very  High. 

the  one  on  energy  demand  is  medium,  then  the  integrated  impact  on  the  calculated  site  scores  that  describe  the  general  risk  level. 

will  be  medium”;  83.60%  selected  “if  the  impact  on  energy  supply 
is  low  and  the  one  on  energy  demand  is  medium-to-high,  then  the 
integrated  impact  will  be  medium-to-high”;  70.50%  selected  “if  the 
impact  on  energy  supply  is  low  and  the  one  on  energy  demand  is 
high,  then  the  integrated  impact  will  be  high”  (Table  3). 

Since  the  integrated  impact  level  can  be  categorized  into  “low”, 
“low-to-medium”,  “medium”,  “medium-high”,  “high”  and  “very- 
high”,  the  corresponding  membership  functions  of  these  fuzzy 
events  can  then  be  established  according  to  Hwang  and  Chen  [62] 
(Fig.  4).  The  range  of  the  integrated  impact  levels  is  subjectively 
assigned  with  a  range  of  [0,1 00].  Thus,  they  will  have  single  numeri¬ 
cal  site  scores  after  defuzzification.  These  numerical  values  have  no 
direct  relationship  with  the  values  of  the  input  impact  factors  (e.g., 
climate  change  impact  levels  of  both  energy  supply  and  demand 
sides).  However,  after  establishing  the  fuzzy  sets  of  the  integrated 
impact  level,  a  numerical  site  score  can  be  obtained  from  Fig.  4 
through  the  fuzzy  “AND”  or  fuzzy  “OR”  operations  based  on  the 
two-side  impact  levels.  This  will  be  illustrated  in  the  following  case 
study.  The  site  management  decisions  can  then  be  made  based 


Table  4  lists  the  relationship  between  site  scores  and  suggested 
management  actions. 

2.3.  Adaptation  planning  under  uncertainty 

The  results  of  integrated  climate  change  impact  analysis  can 
only  provide  information  related  to  “how  serious  will  climate 
change  affect  the  EMS?”  to  decision  makers.  Decision  makers  prob¬ 
ably  need  to  know  “how  and  when  to  take  adaptation  actions 
toward  climate  change?”  Therefore,  it  is  desired  to  couple  the 
impact  analysis  results  with  adaptation  planning  efforts.  Also, 
uncertainties  associated  with  parameter/coefficients  of  adaptation 
planning  need  to  be  dealt  with.  Thus,  an  inexact  optimization 
approach  will  be  adopted.  In  this  research,  an  inexact  energy  model 
(IEM)  will  be  developed  for  supporting  long-term  energy  manage¬ 
ment  systems  planning  under  the  business  as  usual  (BAU)  case, 
as  well  as  adaptation  planning  under  climate  change.  A  number 
of  issues  concerning  cost-effective  allocation  of  energy  resources, 
technologies,  activities  and  investments  can,  thus,  be  addressed. 
It  is  expected  that  the  results  could  provide  useful  information 
and  decision  support  for  energy  management  systems  planning 
within  a  provincial  context.  IEM  is  a  large-scale  dynamic  pro¬ 
gramming  model  for  energy  systems  planning  under  uncertainty 
and  complexity.  This  model  is  subject  to  a  number  of  constraints, 
such  as  end-user  demand,  economic  cost  and  energy  availability.  It 


Table  4 

General  adaptation  strategies  under  coupled  climate  change  impacts. 


Fig.  4.  Membership  functions  of  integrated  climate  change  impact  levels  on  energy 
management  systems. 


Calculated  score 

Adaptation  planning  measures 

90-100 

Immediate  actions  need  to  be  taken  for  mitigating  climate 
change  impacts  on  energy  management  systems 

70-90 

Take  full  actions  for  mitigating  climate  change  impacts  on 
energy  management  systems 

50-70 

Restrict  energy  export  and  consumption 

30-50 

Take  interim  control  measures 

10-30 

Effects  of  climate  change  on  energy  management  systems 
should  be  monitored 
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can  be  used  to  describe  detailed  characteristics  of  energy- related 
technologies  and  evaluate  economic  and  environmental  effects  of 
relevant  policies  and  regulations.  Such  an  energy  management 
system  can  be  characterized  as  a  network/chain  of  energy  flows, 
starting  from  energy-supply  options  (such  as  energy  production, 
export,  and  import)  and  ending  in  end-use  sectors  (such  as  resi¬ 
dential,  commercial,  industrial,  and  transportation  sectors).  It  can 
also  be  customized  into  several  sub-sectors  based  on  site-specific 
features  when  it  is  applied  to  real-world  cases  such  as  in  Province 
of  Manitoba. 

With  consideration  of  local  energy  management  policies,  the 
objective  function  is  to  minimize  the  total  costs  associated  with 
energy  services,  activities  and  investments  in  the  study  province 
over  the  entire  planning  horizon.  This  objective  function  is  a  com¬ 
bination  of  costs  associated  with  a  variety  of  processes  that  move 
energy  from  production  and  conversion  to  consumption.  The  cost 
of  each  process  is  further  divided  into  sub-costs,  such  as  costs 
of  energy  resource  supply,  energy  process/conversion,  technology 
investment,  and  environmental  management.  Decision  variables 
in  IEM  represent  key  discrete  points  of  energy  activities.  Elec¬ 
tricity  generated  from  different  facilities,  production  capacities  of 
energy  generation  facilities,  volumes  of  oil  imported,  and  gaso¬ 
line  consumed  by  private  automobiles  can  be  examples  of  discrete 
variables.  The  constraints  include  the  relationships  among  decision 
variables,  energy  production  efficiencies,  capital  costs,  investment 
costs,  and  energy  availabilities.  Also,  in  order  to  address  uncertain¬ 
ties  associated  with  many  factors/processes  within  an  EMS  and  can 
be  expressed  as  interval  numbers,  an  inexact  programming  method 
is  adopted.  The  model  is  based  on  the  interval-parameter  linear 
programming  (ILP)  model,  in  which  the  uncertain  parameters  are 
expressed  as  intervals  without  any  distributional  information  that 
is  always  required  in  fuzzy  and  stochastic  programming.  The  ILP 
allows  the  interval  information  to  be  directly  communicated  into 
the  optimization  process  and  resulting  solution  [69].  According  to 
Huang  et  al.  [70-72],  an  ILP  model  can  be  written  as  follows: 


±v± 


Min/  =  C±X 
subject  to: 
A±X±  <  B± 


X*  >  0 


(la) 

(lb) 

(lc) 


where  A^IR*}"^",  B^IR*}"^1,  C±e{R±}lxn,  X±e{R±}nxl, 
and  R±  denotes  a  set  of  interval  numbers. 

Thus,  IEM  can  be  formulated  as  follows: 


I  T 


Objective  function  :  Min/*  = 


C; 


l,t 


(2a) 


1=1  t= 1 


( 1 )  The  total  cost  of  energy  consumption  in  the  residential  sub¬ 
sector 


X/u  =  CU  +  Ch2  +  Ct  3  +  Ct,  4  +  Cf,l 


(2b) 


t=l 


(la)  Sub-cost  of  natural  gas  consumption  in  the  residential  sub¬ 
sector 


Cf,l  -  yXu^l.U  +  XU,4,t  +  Xl±l,5,t  +Xf,l,6,t) 


(2c) 


t= l 


(lb)  Sub-cost  of  fuel  oil  consumption  in  the  residential  sub¬ 
sector 


C±  ~~  +  Xf,2,3,t  +  Xf,2,4 


(2d) 


(lc)  Sub-cost  of  liquefied  petroleum  gas  (LPG)  consumption  in 
the  residential  sub-sector 

T 

Cf,3  =  ^^ff,3,t(Xf,3,l,t  +  Xf,3,3,t  +  Xf,3,4,t  +  Xf,3,5,t  +  Xf,3,6,f)(2e) 
t=l 

(ld)  Sub-cost  of  heating  oil  consumption  in  the  residential  sub¬ 
sector 


Cf,4  ~~  ^^fM,t(Xf,4,l,t  +  Xf,4,4,t) 


(2f) 


t= l 


( 1  e)  Sub-cost  of  renewable  energy  resources  consumption  in  the 
residential  sub-sector 


c±  ~~  y^Rtt5,t(Xt5,U  +Xf,5,2,t  +  Xf,5,4,t  +Xf,5,6,t) 


(2g) 


t= l 


(If)  Sub-cost  of  electricity  consumption  in  the  residential  sub¬ 
sector 

T 

Ct,6  =  6,  t(Xf,  6,1  ,t  +Xf,6,3,t  +Xf,6,4,t  +Xf,6,5,t  +  Xf,6,6,tX2h) 

t=l 

(2)  The  total  cost  of  energy  consumption  in  the  commercial  sub¬ 
sector 


+  C±2  +  c2,3  +  C2,4  +  C2,5  +  C2, 


6 


(2i) 


t= 1 


(2a)  Sub-cost  of  natural  gas  consumption  in  the  commercial  sub¬ 
sector 


C2,l  ~~  ^^f2,l,t(X2,l,l,t  +  X2=,l,3,t  +XJ,1,6 


(2j) 


t=  1 


(2b)  Sub-cost  of  fuel  oil  consumption  in  the  commercial  sub¬ 
sector 


C2,2  =  y^f2,2,t(X2,2,l,t  +X^,2,3,t  +X^,2,4,t  +  Xf,2,6,f) 


(2k) 


t=  1 


(2c)  Sub-cost  of  LPG  consumption  in  the  commercial  sub-sector 

T 

C±  =  5^^f2,3,^XJ,3,l,t  +XJ,3,3,t  +XJ,3,5,t  +  Xf,3,6  ^1) 


t= 1 


(2d)  Sub-cost  of  heating  oil  consumption  in  the  commercial  sub¬ 
sector 


C2,4  ~~  ^^f2,4,t(X2,4,l,t  +  X2,4,3,f) 


(2m) 


t= l 


(2e)  Sub-cost  of  renewable  energy  resources  consumption  in 
commercial  sub-sector 


C2,5  _  ^^f2,5,t(X^,5,l,t  +X^,5,3,t  +  X^5,6,f) 


(2n) 


t=  1 


(2f)  Sub-cost  of  electricity  consumption  in  commercial  sub¬ 
sector 

T 

C±  =  ^^f2,6,f(X2,6,l,t +X2,6,2,t +X2,6,3,t +X2,6,4,t +X2,6,5,t 


t= 1 


+  XZ6,6,t) 


(20) 


t= 1 
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(3)  The  total  cost  of  energy  consumption  in  the  transportation 
sub-sector 


~~  C3,l  +  C3,2  +  C3,3  +  C3,4  +  C3,5  +  C3,6  +  C: 


(2p) 


(3a)  Sub-cost  of  natural  gas  consumption  in  the  transportation 
sub-sector 

T 


+  X3,l,5,t  +X3,l,6,t) 


(2q) 


(3b)  Sub-cost  of  fuel  oil  consumption  in  the  transportation  sub¬ 
sector 


C3  2  =  ^  ^3  2  t^X3  2  5  t  +  X3=2  6  t  +  X3=2  7  t ) 


(2r) 


(3c)  Sub-cost  of  LPG  consumption  in  the  transportation  sub¬ 
sector 


C3,3  ~~  YMsMs'U  +XT3,2,t  +X^,3,3,t  +  X3=,3,5,f) 


(2s) 


(3d)  Sub-cost  of  diesel  consumption  in  the  transportation  sub¬ 
sector 

T 

C3,4  =  y^^f3,4,t(X3,4,l,t  +X3,4,2,t  +X3,4,3,t  +  X^4,5,t  +X^,4,6,t^2t^ 
t=  1 

(3e)  Sub-cost  of  jet  fuel  consumption  in  the  transportation  sub¬ 
sector 


C±  “  yZ^3,5,tX3,5,4,t 


(2u) 


(3f)  Sub-cost  of  gasoline  consumption  in  the  transportation  sub¬ 
sector 

T 

CT6  =  y^^f3,6,t(X3,6,l,t  +X^,6,2,t  +  X3=,6,3,t  +X^,6,5,t  +  X3,6,6,t^2v) 
t=  1 

(3g)  Sub-cost  of  electricity  consumption  in  the  transportation 
sub-sector 


(4c)  Sub-cost  of  LPG  consumption  in  the  industrial  sub-sector 


M  T 


C±  — 
l4,3  ~~ 


D/-±  y± 

m4,3,tA4,3,m,t 


(2aa) 


m  t=  1 


(4d)  Sub-cost  of  coal  consumption  in  the  industrial  sub-sector 

T 

C4,4  =  y^f4,4,f(X4,4,l,t  +X4,4,3,t  +  X4,4,5,t  +  X4,4,6,f)  (2ab) 

t=  1 

(4e)  Sub-cost  of  renewable  energy  resources  consumption  in  the 
industrial  sub-sector 


M  T 


C±  — 

L4,5  “ 


D/-±  y± 

AL4,5,tA4,5,m,t 


(2ac) 


m  t=  1 


(4f)  Sub-cost  of  electricity  consumption  in  the  industrial  sub¬ 
sector 


M  T 


C±  — 
l4,6  ~~ 


D/-±  y± 

A'L4,6t'A4,6,m,t 


(2ad) 


m  t= 1 


(5)  The  total  cost  of  power  generation  facilities 


N  T 


J2C±  =  [(AF0Mn,t  +  AVOM£tY*  )  +  CAP±t  INV^Z^; 


n= 1  t= 1 


(2ae) 


subject  to: 

(1 )  Balance  of  energy  carriers 
(la)  Balance  of  natural  gas 


2  5 


,m,t  + 


2=1  m= 1 


m= 1 


m= 1 


Ele± 

Conv? 


^  <^tIMP-t,  Vt 


(2af) 


(lb)  Balance  of  fuel  oil 


kf,2,l,t  +Xf,2,3,t  +Xf,2,4,t  +X^,2,l,t  +X^,2,3,t  +X^,2,4,t  +X^,2,6,t 


+  X3,2,5,t  +X3,2,6,t  +X^,2,7,t  +  y^X4,2,m,t  -  ^2,t^P2,r 


m~  1 


(2ag) 


C3.7  ~~  +  X3J,3,t  +X3J,4,t  +X3,7,5,t  +  X3,7,6,t)(2w)  Ra*anCe  °f  LPG 


(4)  The  total  cost  of  energy  consumption  in  the  industrial  sub¬ 
sector 


t  =  Ct  1  +  C4,2  +  CJ,3  +  Ct,4  +  Ctl 


(2x) 


(4a)  Sub-cost  of  natural  gas  consumption  in  the  industrial  sub¬ 
sector 


M  T 


C±  — 
L4,l  - 


D/-±  y± 

m4,l,tA4,l,m,t 


(2y) 


m  t= 1 


(4b)  Sub-cost  of  fuel  oil  consumption  in  the  industrial  sub-sector 


c±  — 

l4,2  ~~ 


D/-±  y± 

AL4,2,tA4,2,m,t 


(2z) 


Xf,3,l,t  +Xf,3,3,t  +Xf,3,4,t  +Xf,3,5,t  +Xf,3,6,t  X2=,3,l,t  +XJ,3,3,t 

+  X^,3,5,t  +X^,3,6,t  +X^,3,l,t  +X^,3,2,t  +X3=,3,3,t  +X^,3,5,t 


M  T 


(2ah) 


m  t=l 


(ld)  Balance  of  heating  oil 

xi,4,i,t+XiMt+^U+^4,3,t  ^  Vf  (2ai) 

(le)  Balance  of  renewable  energy  resources 

Xf,5,l,t  +Xf,5,2,t  +Xf,5,4,t  +Xf,5,6,t  +X2,5,\,t  +XJ,5,3  ,t  +  XJ,5,6  ,t 


M  T 


Ks,tKs,m,t  *  DREN?  •  Vt 


(2aj) 


m  t= 1 


m  t=  1 
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(If)  Balance  of  electricity 


(2h)  Freight  truck  energy  demand 


2V  j 

EleNGA*  +  EleFO*  +  EleCOALt'  +  EleRENt'  >  V\y*t),  Vt(2al<)  ^Con vfj  ^j  2  [  >  DFR'f,  Vi  =  3,  tj  =  1, 3,  4,  6and7  (2aw) 


(lg)  Balance  of  diesel 

X  ^  3  |  f  “F  ^3~  3  2  f  T  3  3  f  "I"  3  5  f  “F  ^3”  3  g  f 

Y± 

+  Vt 

ConvJt 

(lh)  Balance  of  jet  fuel 

^4^i1MP6,t-  Vf 

(li)  Balance  of  gasoline 

^3,6,l,t  +  ^3,6,2,t  +  ^3,6,3,t  +  ^3,6,5,t  +  ^3,6,6,t  - 


(lj)  Balance  of  coal 


Y± 

r3  ,t 


(2al) 


(2am) 


(2an) 


Xf,4,U  +  X4±4,3,t  +  X4,4,5,t  +^4,6,t  +  Vt 


(2ao) 


(2i)  Bus  energy  demand 


^Conv*  fX*  3  f  >  DBS *  ,  Vi  =  3,  t,  j  =  1 , 3, 4,  6  and  7  (2ax) 


(2j)  Air  plane  demand 


EConvU,tX5,4,t  *  D/VR±,  Vi  =  3,  t,  j  = 


(2k)  Recreational  energy  demand 


(21)  Railway  energy  demand 


(2ay) 


EConv*  X*  ,  f  >  DRC*  Vi  =  3,  t,j  =  1, 3, 4,  6  and  7  (2az) 

i,2,r  1 J  >  J ,  r 


^7conv*.  fX*.  g  f  >  DRL-'f,  Vi  =  3,  tj  =1,2,  4,  6,  7  for  railway 


(2ba) 


(2)  End-user  energy  demand  constraints 
(2a)  Space  heating  demands 

J 

ECon*  X*.  ,  ,  >  DSH*  Vi,  ti  =  land2 

l,J,t  lj,l,t  -  l,t’ 

2=1 

(2b)  Space  cooling  demands 
J 

ECon*.  X1.  ,  ,  >  DSC'  ,  Vi,  t,  i  =  1  and 2 
2=1 

(2c)  Water  heating  demands 
J 

V Con*.  .X*  ,  .  >  DWH*  ,  Vi,  t,  i  =  1  and  2 

/  -J  hj  h  j  9  ^  9  ^  *  9  ^ 

2=1 

(2d)  Lighting  demands 
J 

ECon1  X*.  .  „  >  DLG '  ,  Vi,  t,  i  =  1  and  2 

i,j,t  lj, 4, t  -  i,V  ’  ’ 

2=1 

(2e)  Appliance  energy  demands 
J 

V Con*  .X*  ,  .  >  DAP*  ,  Vi,  t,  i  =  1  and  2 
2=1 

(2f)  Commercial  ventilation  demand 
J 

ECon*.  X*  .  „  >  DVP*  ,  Vi,  t  i  =  2 

1,2, t  1,2, 4, t  -  i,V  ’ 


(2ap) 


(2aq) 


(2ar) 


(2as) 


(2at) 


(2au) 


(2m)  Equivalent  lubricant  energy  demand 


EConvf.  Xf.  7  >  DLB^,,  Vi  =  3,  t,  j  =  2  for  lubricant  (2bb) 

*  9  J  9  ^  ^  9  '  9  ^  l  y  L 


(2n)  Total  electricity  demand 


2  5 


^i,6,m,t  ^  7^3,6,m,f  ^  y^4,6,m,t 


i=l  m=l 


171  =  1 


171  =  1 


+  EELEf  <t)*t^Y*t,  Vt 


(2o)  Agricultural  energy  demand 


EConV«,^U  ^  DAGCM’  V/  =  4,  t,  j  =  1, 2,  . . . ,  6 


(2p)  Forest  energy  demand 


EC0I1VmA/.:u  *  DFR05-  Vi'  =  4,  t,  j  =  1 , 2,  3,  5,  6 


(2q)  Mining  energy  demand 


EConVU.,XU.  u  *  DMINtr  Vi  =  4,  t,j  =  1, 2, . . . ,  6 


(2r)  Construction  energy  demand 


(2bc) 


(2bd) 


(2be) 


(2bf) 


ZConVU.,XU,u  ^  DC0NM’  Vi  =  4,  t,  j  =  1, 2,  3,  5,  6  (2bg) 


(2g)  Light  duty  vehicle  energy  demand 


(2s)  Pulp  and  paper  energy  demand 


Jconv*  fX*  ,  f  >  DLD*  ,  Vi  =  3,  t,j  =  1,3,4,  6and7  (2av)  J^Conv*  fX*  5  f  >  DPP*  ,  Vi  =  4,  t,  j  =  1, 2, . . . ,  6 


(2bh) 
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(2t)  Manufacture  energy  demand 


3.  Method  of  solution 


J 


EConvu ,rXu ,e ,t  ^  DMANu-  Vi  =  4,  t ,j  =  1 , 2, . . . ,  6  (2bi) 


(2u)  Electricity  export  targets 
EELE*  >  TARGET_EX*,  Vt 

(3)  Electricity  generation  facility  capacity  constraint 
Yn,t  5  EXCAP, tf  +  CAP^Z^ 

Lower  Y*t  <  Y*t  <  Upper Y*t,  Vn,  t 

(4)  GHG  emission  constraints 


/  J 


M 


(4a)EEEmFtu,fE^,,f  2  (EmC0^ 

i=l  j=l 


/  J 


m=l 


M 


(4b)EEEmFt2,,,fExu.m.f  *  (£mCH^ 


i=l  J=1 


/  J 


m~  1 


M 


(4c)EEEmFt3,'J,fEE,.,f  ^  (EmN2°)f 


i=l  J=1 


m=l 


(2bj) 


(2bk) 

(2bl) 


(2bm) 


(2bn) 


(2bo) 


where  t  =  1 ,2,. .  .,5;  Ci  i  =the  total  cost  of  natural  gas  (NGA)  con¬ 
sumption  in  residential  sub-sector;  c12=the  total  cost  of  fuel  oil 
(FO)  consumption  in  residential  sub-sector;  c\ >3  =  the  total  cost  of 
liquefied  petroleum  gases  (LPG)  consumption  in  residential  sub¬ 
sector;  ci>4  =  the  total  cost  of  heating  oil  (HO)  consumption  in 
residential  sub-sector;  ci>5=the  total  cost  of  renewable  energy 
(RE)  resources  consumption  in  residential  sub-sector;  c16  =  the 
total  cost  of  electricity  (ELE)  consumption  in  residential  sub¬ 
sector;  c2ii=the  total  cost  of  NGA  consumption  in  commercial 
sub-sector;  c2>2  =the  total  cost  of  FO  consumption  in  commercial 
sub-sector;  c2>3  =the  total  cost  of  LPG  consumption  in  commercial 
sub-sector;  c2>4  =  the  total  cost  of  HO  consumption  in  commercial 
sub-sector;  c2  5  =the  total  cost  of  RE  consumption  in  commercial 
sub-sector;  c2  6  =  the  total  cost  of  ELE  consumption  in  commercial 
sub-sector;  c31=the  total  cost  of  NGA  consumption  in  trans¬ 
portation  sub-sector;  c3  2  =  the  total  cost  of  FO  consumption  in 
transportation  sub-sector;  c3>3  =  the  total  cost  of  LPG  consump¬ 
tion  in  transportation  sub-sector;  c34  =  the  total  cost  of  diesel 
(DSL)  consumption  in  transportation  sub-sector;  c35=the  total 
cost  of  jet  fuel  (JF)  consumption  in  transportation  sub-sector; 
c3  6  =  the  total  cost  of  gasoline  (GSL)  consumption  in  transportation 
sub-sector;  c37  =  the  total  cost  of  ELE  consumption  in  trans¬ 
portation  sub-sector;  c44=the  total  cost  of  NGA  consumption 
in  industrial  sub-sector;  c42=the  total  cost  of  FO  consumption 
in  industrial  sub-sector;  c43=the  total  cost  of  LPG  consumption 
in  industrial  sector;  c44  =  the  total  cost  of  coal  (CL)  consump¬ 
tion  in  industrial  sector;  c45=  total  cost  of  RE  consumption  in 
industrial  sector;  c4>6  =  total  cost  of  ELE  consumption  in  indus¬ 
trial  sector;  n  =  electricity  generation  facilities  (n  =  1  =  hydropower; 
n  =  2  =  coal;  n  =  3  =  NGA;  n  =  3  =  fossil  fuel);  rjf  t  =  coefficient  of  loss; 

rjg  t  =  line  loss  of  electricity  transportation. 

Obviously,  the  above  IEM  model  can  effectively  reflect:  (i)  inter¬ 
actions  between  a  number  of  sub-sectors,  technologies  and  energy 
resources  among  various  energy  management  systems,  (ii)  uncer¬ 
tain  information  associated  with  input  parameters  that  can  be 
expressed  as  intervals,  and  (iii)  the  provincial  baseline  of  energy 
management  over  the  planning  horizon.  Thus,  a  series  of  decision 
related  to  “how  and  when  to  take  adaptation  actions  toward  climate 
change”  can  be  obtained. 


The  method  of  solution  is  based  on  (a)  the  results  of  FIIM  for 
climate  change  impact  analysis,  and  (b)  the  interactive  solution 
algorithm  of  ILP.  Such  an  interactive  solution  algorithm  has  been 
developed  to  solve  model  (1 )  through  analyzing  the  detailed  inter¬ 
relationships  between  parameters  and  variables  and  between  the 
objective  function  and  constraints.  According  to  the  algorithms, 
the  solution  for  interval  linear  programming  (ILP)  model  can  be 
obtained  through  a  two-step  method,  where  a  submodel  corre¬ 
sponding  to  f~  (when  the  objective  function  is  to  be  minimized) 
is  first  formulated  and  solved,  and  then  the  relevant  submodel  cor¬ 
responding  to  f  can  be  formulated  based  on  the  solution  of  the  first 
submodel.  In  detail,  the  first  submodel  can  be  formulated  as  follows 
(assume  that  bf  >  0,  and/*  >0): 


n 


Min/'  =  EC7*T  +  E  CiX! 

j= 1  j=k  i+l 

subject  to: 


(3al) 


n 


Ejai)l+Sign(ai;  )x;  +  E  '“'J'1  SiSn(ay)X/  ^  bt ’ 

j= l  j=k  i+l 


Vi  (3a2) 


xf  >  0,  Vj  (3a3) 

where  x*,  j  =  1,2,. .  .,/<i  are  interval  variables  with  positive  coef¬ 
ficients  in  the  objective  function  and  x*,  /  =  /<i  +  l,/<i  +2,. .  .,n  are 
interval  variables  with  negative  coefficients  in  the  objective  func¬ 
tion.  Thus,  solutions  of  xj  {j=  1,2,. .  .,k\)  and  ^+opt  (j  =  /<i  + 1 , 
k\  +2,. .  .,n)  can  be  obtained  through  solving  submodel  (3a).  Then 
the  submodel  corresponding  to  f+  can  be  formulated  as  follows 
(assume  that  bf  >  0,  and/*  >  0): 


n 


Min/*  =  c*x*  +  c*x. 

/  j  J  J  A _ J  J  J 

J=l  j=k  i+l 

subject  to: 


(3bi) 


n 


^Ja0|  Sign (aiJ)xf  +  ^  |a0|+Sign(at)x.  <bf, 
j= l  j=k  i+l 


Vi 


(3b2) 


xf  >  0,  Vj  (3b3) 

xf  >  xj opt,  j  =  1,2,  (3b4) 

xj  <  xf  opt,  j  =  ki  +1,  k2+2,...,n  (3b5) 

Solutions  ofx/opt  0  =  1.2 . fc, )  and  xT  opt  (j  =  ki  +1/1  +2 . n) 

can  be  obtained  through  solving  submodel  (3b).  Thus,  we  can  have 
the  final  solution  of ffpt  =  and  xf  opt  =  [xj  opt,  xf  opt]. 

Thus,  the  integrated  modeling  system  (IEMS)  can  be  developed 
through  linking  the  fuzzy-interval  inference  method  (FIIM)  and  the 
inexact  energy  model  (IEM).  The  general  solution  procedure  for  the 
proposed  IEMS  is  presented  in  Fig.  5  and  described  as  follows: 


-  Step  1:  Identify  climate  change  impacts  on  energy  supply  based 
on  the  obtained  fuzzy  sets; 

-  Step  2:  Identify  climate  change  impacts  on  end-user  demand 
based  on  the  obtained  fuzzy  sets; 

-  Step  3:  Identify  multi-layer,  multi-sector  and  multi-period 
impacts  of  climate  change  on  the  EMS  through  the  integration 
method; 

-  Step  4:  Based  on  interval  analysis  and  fuzzy  sets  theory,  identify 
the  integrated  impact  levels  of  climate  change  on  the  EMS; 
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Climate  change  impacts 


Decomposing 


Energy  supply 

Energy  End-user 

<Zfy$ 

43-  Coupling 

Multi-layer, 

multi-sector  and 

multi-period  impacts 

Interval  analysis  Fuzzy  sets  theory 


Integrated  impact  analysis 


43- 


Interval  solutions  (decision  alternatives) 


Fig.  5.  Solution  procedure. 


-  Step  5:  Based  on  the  integrated  impact  results  and  the  values  in 
Table  4,  judge  whether  adaptation  actions  would  be  needed  for 
the  EMS  under  climate  change; 

-  Step  6:  If  no  adaptation  actions  would  be  needed,  formulate  the 
inexact  model  (IEM)  under  the  BAU  case; 

-  Step  7:  Transform  the  IEM  model  into  two  submodels,  where  the 
lower  bound  of fi=  is  first  desired  since  the  objective  is  to  minimize 

J*; 

-  Step  8:  Formulate  the  f~  submodel,  including  the  f~  objective 
function  and  the  related  constraints; 

-  Step  9:  Solve  the/-  submodel  according  to  conventional  LP  (linear 
programming)  solution  methods,  and  obtain  the  solution  for f~t 

andxr  • 

I  opt’ 

-  Step  10:  Formulate  the/+  submodel,  including  the  objective  func¬ 
tion  and  constraints  corresponding  to/+,  as  well  as  the  second  set 
of  constraints  regulating  the  bounds  of  decision  variables; 

-  Step  11:  Solve  the/+  submodel  according  to  conventional  LP  (lin¬ 
ear  programming)  methods  and  obtain  the  solution  for/0+t  and 

y  H-  • 

j  opt’ 

-  Step  12:  The  optimal  results  for  the  IEMS  model  can,  thus,  be 
obtained  as  =  Lfopt,/0+pt]  and  xf  opt  =  [xj  opt,  x+  opt]; 

-  Step  13:  If  adaptation  actions  need  to  be  undertaken  in  respond¬ 
ing  to  the  integrated  climate  change  impacts,  formulate  the 
corresponding  model  based  on  modeling  structure  formulated 
under  the  BAU  case  and  the  impact  analysis  results  (i.e.,  the  vari¬ 
ations  on  energy  supplies  and  demands)  obtained  from  Steps  1  to 
5; 

-  Step  14:  Repeat  Steps  7  to  12; 

-  Step  15:  Make  comparisons  between  results  under  the  BAU  case 
and  climate  change. 


4.  Conclusions 

In  this  research,  a  large-scale  integrated  modeling  system  (IMS) 
was  developed  for  supporting  adaptation  planning  under  impacts 


of  climate  change  from  an  entire  EMS  point  of  view.  A  number 
of  methodologies  were  seamlessly  integrated  into  IMS,  includ¬ 
ing  the  fuzzy-interval  inference  method  (FIIM),  inexact  energy 
model  (IEM),  and  uncertainty  analysis.  Compared  to  conventional 
studies  in  the  similar  areas,  such  a  system  could  (i)  incorporate 
multiple  technologies,  energy  resources,  and  sub-sectors,  and  cli¬ 
mate  change  impacts  within  a  general  modeling  framework,  (ii) 
address  interactions  of  climate  change  impacts  on  multiple  energy 
sub-sectors  and  resources  within  an  EMS,  (iii)  identify  optimal 
adaptation  strategies  of  an  EMS  to  climate  change  impact  through 
a  two-step  procedure,  (iv)  deal  with  multiple  levels  of  uncertain 
information  associated  with  processes  of  climate  change  impact 
analysis  and  adaptation  planning,  and  (v)  seamlessly  combine 
climate  change  impact  analysis  results  with  inexact  adaptation 
planning.  Thus,  decision  makers  could  have  a  comprehensive  view 
on  the  EMS  as  well  as  the  corresponding  adaptation  schemes  under 
climate  change,  greatly  improving  the  robustness  and  complete¬ 
ness  of  the  decision-making  processes.  The  developed  method  can 
be  applied  to  long-term  adaptation  planning  in  energy  manage¬ 
ment  systems,  particularly  those  that  are  highly  dependent  on 
renewable  energy  sources.  The  solutions  can  generate  decision 
alternatives  and  thus  help  decision  makers  identify  desired  policies 
under  climate  change.  In  the  companion  paper,  application  of  the 
developed  approach  will  be  conducted  in  the  Province  of  Manitoba, 
Canada. 
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